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Abstract 

The matrix elements of the multi-channel Jost matrices are written in such a way 
that their dependencies on all possible odd powers of channel momenta are factorized 
explicitly. As a result the branching of the Riemann energy surface at all the channel 
thresholds is represented in them via exact analytic expressions. The remaining single- 
valued functions of the energy are expanded in the Taylor series near an arbitrary point 
on the real axis. Using the thus obtained Jost matrices, the S'-matrix is constructed and 
then the scattering cross section is calculated, which therefore depends on the Taylor 
expansion coefficients. These coefficients are considered as the adjustable parameters 
that are optimized to fit a given set of experimental data. After finding the coefficients, 
the resonances are located as zeros of the Jost matrix determinant at complex energies. 
Within this approach the S'-matrix has proper analytic structure. This enables us not 
only to locate multi-channel resonances but also to reproduce their partial widths as 
well as the scattering cross section in the channels for which the data are not available. 
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1 Introduction 



Quantum resonances play an important role in theoretical description and intuitive understand- 
ing of physical processes taking place in the microscopic world of molecules, atoms, nuclei, and 
various nano-systems. This is why locating the corresponding spectral points £ = — |r, 
i.e. determining the resonance parameters, Er (resonance energy) and V (resonance width), 
is an indispensable part of any theoretical modelling as well as of an analysis of scattering data. 

From the theoretical point of view, the resonances, being complex-valued eigen-energies of 
the Hamiltonian, are not easy to locate. For each angular momentum, they appear as regular 
strings on the so called non-physical sheet of the complex energy surface (see, for example, 
Ref. [Ij). The influence of these resonances on the corresponding partial-wave cross section 
depends not only on their distribution over the energy surface but also on the corresponding 
5*- matrix residues [2i-i5j. 

The analysis of scattering data, i.e. phenomenological extraction of resonance parameters, 
is perhaps even more difficult than their theoretical determination. Indeed, a measured cross 
section is the sum of all partial-wave cross sections, where the effects of all possible reso- 
nances add together. Therefore the sharp structures corresponding to individual resonances, 
may overlap and smear out their discrete character. Such an overlap very often happens even 
with resonances belonging to the same partial wave. 

In a phenomenological search for resonances, the first step usually consists in performing the 
so called phase-shift analysis [HEllZ], i e. in obtaining numerical values of the partial-wave 
phase shifts or the corresponding reaction amplitudes, for a set of collision energies. Then 
the amplitudes known on the real axis, are analytically continued onto the complex-energy 
Riemann surface, where the resonances are identified as their singularities. Various methods 
differ in the way such an analytic continuation is done pLj. 

The basic idea underlying the majority of existing approaches of this type, consists in fitting 
the experimental points with curves obtained from a phenomenological scattering amplitude, 
where possible resonance poles are artificially embedded by hand and their complex energies 
serve as the fitting parameters. The most simple of such approaches is the Breit-Wigner 
parametrization [8J. Another version of the same approach is based on the Fano parametriza- 
tion [5], where the amplitude is split in the resonant and background parts, which makes it 
more realistic and allows one to treat more complicated "zigzags" of the cross section. In nu- 
clear and particle physics, the resonances are usually introduced into the fitting procedure via a 
model propagator with explicit singularities at complex energies (unitary isobar model) [TUlfTI] . 
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Many authors emphasize (see, for example, Chapter 6 of the book [T] and references therein as 
well as a more recent publication [12J) that for a reliable analytic continuation of the scattering 
amplitude (or S-matrix) it is very important to construct its phenomenological expression pre- 
serving proper analytic structure, i.e. its branching points, various symmetries etc. This can 
easily be achieved for the single-channel problems and with some effort for two-channel ones 
(Dalitz-Tuan representation) [12]. However, when the number of channels is three or more, 
constructing the S-matrix with correct branching at all the thresholds becomes extremely dif- 
ficult within the traditional approaches [131 [14] . 

In the present paper, we resolved this difficulty. For obtaining the S'-matrix, we use the 
structure of the Jost matrix derived in our early publications [l5l[T5|, where all the factors 
responsible for the branching are given analytically for an arbitrary number of channels. This 
enables us to perform the analytic continuation and the search for resonances in a model- 
independent way, i.e. without embedding the resonance poles by hand. 



2 Jost matrices 

Consider an A^-channel collision process that can symbolically be written as a chemical or 
nuclear reaction, 

' A + B , 
A + B ( C + D , 



with two-body systems in both the initial and final states. We assume that the interaction 
forces for these systems are of a short-range type, i.e. they vanish at large distances (r — )• oo) 
faster than the Coulomb force. The wave function describing such a process at the collision 
energy is a column matrix, 

^{E,r)=\^^{E,r) , 

where each of the lines corresponds to a separate channel. The channels are characterized 
not only by the type of the particles [A + B, C + D, etc.) but also by the complete sets of 
the corresponding quantum numbers (such as the threshold energy, angular momentum, spin. 
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etc.). In other words, any two states of the system that differ by at least one quantum number 
are considered as different channels even if the type of the particles and the threshold energies 
are the same. This means that a channel state has a definite value of the angular momentum 
and thus its angular dependence can be factorized as 

r 

The set of radial wave functions u„ obey the coupled system of radial Schrodinger equations, 



+ kl 



Un{E,r) = ^ ^ Vnn'{r)Un'{E,r) 



N 



(1) 



n'=l 



where the coupling is due to the off-diagonal elements of the interaction-potential matrix V^, 
In Eq. ([T]), the channel momentum. 



2/iri 



{E - En) 



(2) 



is determined by the difference between the total energy E and the channel threshold En, as 
well as by the reduced mass fin in the channel n. 

Each of the equations of the set ([T]) is of the second order. In the theory of differential 
equations (see, for example, Ref. [17j) it is shown that such a set has 2N linearly independent 
solutions, i.e. 2N linearly independent columns that form a basis in the solution space, and only 
half of these columns are regular at r = 0. The regular columns can be combined in a square 
matrix ^{E,r), which is called the fundamental matrix of regular solutions. Since a physical 
wave function must be regular, it is a linear combination of the columns of the matrix ^{E, r). 

When the particles move away from each other (r — )• oo), the potential matrix tends to zero 
and thus the right hand side of Eq. ([T]) vanishes. The remaining set of N second order 
(Riccati-Bessel) equations, 



dy, ~1~ k^ 



Un{E,r) 







r — > oo 



(3) 



has 2N linearly independent column-solutions. These 2N columns can be combined in two 
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diagonal square matrices, 



p^(out) 





V 



V 













\ 





\ 




(4) 



(5) 



that involve the Riccati-Hankel functions h^^\kr) and represent the m-coming and owt-going 
spherical waves in all N channels. 



The 2N columns of the matrices (Hj) and ([5]) constitute a basis in the solution space at large 
distances and thus each column of the fundamental matrix r) is their linear combination 
when r — 7- oo. This can be written as 

$(E,r) — ^ iy('"^)(E,r)F('"^)(E) + iy(°"*)(E,r)F(°'^*)(E) , (6) 

r— >oo 

where the combination coefficients are combined in the square matrices F'^™\E) and F^°^^\E). 
They are functions of the energy and are called the Jost matrices. 

The ^-matrix that completely determines all the scattering observables, is expressed via the 
Jost matrices (the details can be found, for example, in the textbook or in the papers 

S{E) = . (7) 

The resonances are those spectral points 

S = Er-^T , Er>0 , r>0 , (8) 

on the Riemann surface of the energy, where 

detF('")(^) = . (9) 
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The energy surface has a square-root branching point at every channel threshold En. This is 
because the Jost matrices depend on the energy E via the channel momenta ((21) and for each 
of them there are two possible choices of the sign in front of the square root. The resonance 
spectral points are located on the so called non-physical sheet of this Riemann surface, i.e. 
such a layer of the surface where all the channel momenta have negative imaginary parts. In 
the numerical calculations, the choice of the sheet is done by an appropriate choice of the 
signs in front of the square roots ([2]). 

3 Analytic structure of the Jost matrices 

in the present paper, our main goal is to find a way of a reliable parametrization of experimental 
data, such that it would allow us to locate the resonance spectral points at complex energies. 
As we described in the previous section, the multi-channel ^-matrix has very complicated en- 
ergy dependence via the channel momenta and therefore is defined on a Riemann surface with 
an intricate connection of many layers. This means that a straightforward parametrization of 
such a matrix using an arbitrarily chosen functional form may give erroneous results. When 
choosing the parametrization form, it is important to take into account as much information 
on the symmetry properties and analytic structure of the ^-matrix as possible. 

First of all, we notice that the S'-matrix is a kind of "ratio", given by Eq. ([7]), of the Jost 
matrices and which are not completely independent of each other but rather are 

somehow related. In other words, the parameters in the "numerator" and "denominator" of 
(I7D should be the same. Indeed, as we found in Ref. [15], for the systems interacting via 
short-range potentials, the Jost matrices have the following structure: 

F^^iE) = l^A^UE)- BUE), (10) 

^^n\E) = ^^,AUE)+ BUE), (11) 

where the energy dependent matrices A{E) and B{E) are the same for both F^'"^ and 
Moreover, in the same Ref. [16J it was established that the matrices A(E) and B{E) are 
single-valued analytic functions of the energy defined on a single one-layer energy plane. In 
other words, the matrices A{E) and B(E) are the same for all the sheets of the Riemann 
surface and all the complications stemming from the branching points are isolated in Eqs. 
( ITOl) and ( ITTI) via the explicit factors depending on the channel momenta. Therefore, by using 
the expressions ( ITOl ) and ( ITII ) in Eq. we guarantee that the "numerator" is properly 
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related to the "denominator" (which means guaranteed unitarity on the real axis) and that all 
the branching points are properly embedded. If we find an adequate parametrization of the 
matrices A(E) and B{E), the resulting ^-matrix will automatically have correct values on all 
the sheets of the Riemann surface. 

Since the matrices A{E) and B{E) are analytic functions of the variable E, they can be 
expanded in the Taylor series, 

oo oo 

A{E) = Y,iE - EoTa^iEo) , B{E) = J^iE ~ EoTb^Eo) , (12) 

n=0 n=0 

near an arbitrary point Eq within the domain of their analyticity. Here the expansion coeffi- 
cients ttn and bn are matrices of the same dimension as A and B. They depend on the choice 
of the point Eq. 

In Ref. [16J it was shown that for a given potential the expansion coefficients a„ and 6„ can 
be obtained as the asymptotic values of the solutions of certain set of differential equations. 
In the present paper we are not going to calculate these coefficients since we assume that 
the potential is not known. Instead, we will use a„ and bn as fitting parameters in order to 
reproduce experimental scattering cross section. As soon as the optimal expansion coefficients 
corresponding to the experimental data are found, we can use them to obtain the matrices 
A{E) and B{E) and through them the Jost matrix F^'''\E) given by Eq. (ITOl) . Then we can 
locate the resonances as the complex roots of Eq. (0). 

Although we are not going to calculate the expansion coefficients, the fact that they are solu- 
tions of certain differential equations derived in Ref. [16J is important. Indeed, that differential 
equations have real boundary conditions and all their coefficients are real if Eq is chosen on 
the real axis. This means that for real Eq the matrices a„ and 6„ do not have imaginary parts, 
i.e. the total number of fitting parameters is 2(M + 1)A^^ where is the matrix dimension 
(number of channels) and M is the highest power in the expansions 

M M 

A{E) ^ Y.^E - EoTa^iEo) , B{E) ^ ^(^ - ^o)"&n(^o) • (13) 

n=0 n=0 

If the central point Eq of the expansion is taken not on the real axis, then the number of the 
parameters we have to fit is doubled. 

The approximate expansions (TT^ can only be accurate within a circle around Eq. In Ref. [T5] 
it is demonstrated how the radius of such a circle increases with the number of terms taken 
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into account. When using this approach to extract the resonance parameters from experimen- 
tal cross section, we therefore should not fit the data within a too wide energy-interval. The 
central point Eq must be placed somewhere in the middle of the interval where it is expected 
to find a resonance. An adequate width of the interval around Eq unfortunately cannot be 
estimated. It is necessary to repeat the analysis with several different widths and make the 
decision on the basis of the stability of the results thus obtained. 



4 Fitting procedure 

Suppose we have sets of experimental data for several channels, m — t- m' — t- n', etc. 



E. 



(mn) 



± s. 



(mn) 



12 3 Ar(™") 

i, U, . . . , -1 V , 



O'm'n' ( Ej 



{m'n') 



± 6 



(m'n') 



3 

etc. 



i = 1, 2, 3, . . . , A^^*"'"') 



Here amn and S^^"^ are the cross section in the channel m ^ n and the corresponding standard 
deviation while the energy interval covered by the points E^[^^^ is around the energy E^ where 
we expect to find a resonance. In order to parametrize the Jost matrix, we construct the 
function 



X 



i=l 



/p(mn)x fit /p(r7^n)^ 



(mn) 



(14) 



+ E 



am'n'[Ej )-a^,^,{E 



-,{m'n') y 



(m'n') 



+ 



etc. 



where the fitting cross section for the channel m n 



7r 



= — (2C + 1) \Snm{E) - 6„ 



(15) 



depends on the expansion coefficients of (llSp via Eqs. (ITU|) . (ITTI) . and ([7]). These coefficients 
therefore serve as the fitting parameters. 



8 



The time-reversal invariance leads to the so called detailed balance theorem which means that 
the S'-matrix is symmetric with respect to the transposition, i.e. Smn = Snm- 'f we simply min- 
imize the x^-function given by Eq. ( 1141) . this symmetry is not guaranteed. If is the number 
of channels (the dimension of the matrices), then the symmetry gives us (N'^ — N)/2 equa- 
tions (the number of the elements above the diagonal of the matrix) relating the variational 
parameters. By solving this set of equations, we can reduce the number of such parameters. 
Although the equations are not simple, this can always be done numerically for any reasonable 
value of A^. 

There is a more simple way of making the S'-matrix symmetric although it requires to vary a 
bit more parameters. A set of the optimal parameters Oq, ai, . . . , om, ^i; • • • ? i>M that give 
a symmetric S'-matrix, can be obtained by minimizing the generalized x^-function 

'^^(Oo, Ol, • • • , Om, ^1, • • • , ^Af) = + ^ \Smn{Ej) — Snm{Ej)\^ , (16) 

m<n,j 

where at all experimental points the differences between the off-diagonal elements are included. 

After finding the optimal parameters, we obtain analytic expressions for the Jost matrices and 
the S'-matrix, valid within a circle around Eq on all the sheets of the Riemann surface. Using 
these expressions, we should be able not only to locate the nearest resonances but also to 
calculate the cross sections in all the other channels for which we do not have experimental 
data. 

5 Examples 

The proposed procedure for parametrizing experimental cross section and thus locating the 
resonances needs to be demonstrated by a couple of simple and clear examples. In such ex- 
amples, we should know the resonance parameters beforehand. This will give us the feeling of 
how accurate the procedure is. 

As the examples, we chose the well-known and well studied one- and two-channel models 
specified by certain potentials (see the next two sections). We use these potentials to artificially 
generate "experimental" data points around the energies where they support resonances, and 
then try to recover these resonances using the suggested parametrization method. Since the 
exact values of the resonance parameters are known, this shows us how reliable the proposed 
method is. 
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5.1 Single- channel model 



The simple potential barrier given by 

V{r) = 7.5r^e-' , (17) 

is very often used as a testing ground for new theoretical methods [21] ■ In this model, the 
units are such that h^/ fi = 1 and thus the energy as well as the distances are dimensionless. 
It has a rich spectrum of resonances. The first two of them with £ = are [25]: 

^exact ^ 3.426390 - ^0.025549 , (18) 
^exact _ 4.834807-^2.235753. (19) 

Let us assume that we are given experimental S-wave scattering cross section for a single- 
channel system in the energy interval 3 < E < 4.5 as is shown in Fig. [H Actually, these 
15 points are generated using the potential (ITTf) . along the corresponding exact cross section 
shown in Fig. |2](thin curve), which has a sharp zigzag near the first very narrow resonance (IT5I) . 

In order to fit the data, we use the approximate expressions (IT3i) with Eq = 3.4 and M = 5. 
Being substituted into Eqs. (ITU!) and (ITTi) . they give us the Jost functions, from which we 
obtain the ^-matrix (JT]) and finally the cross section (IT5l) . As we mentioned before, for a real 
£"0 the coefficients (fitting parameters) a„ and 6„ are also real. This means that we have to 
adjust 12 parameters that minimize the x^-function of the type (IT4|) where only single elastic 
channel is taken into account. 

As the minimization tool, we used well known program "MINUIT" from the CERN library [221 
UA\ . In order to avoid the situation when we stuck in a local minimum, we repeated the 
minimization procedure several hundreds of times with randomly chosen initial values of the 
parameters and took the best findings. The best minimum we found was with = 1-5 x 10~^. 
The result of this fitting is shown in Fig. [2](thick curve). There is no visible difference between 
the exact (thin curve) and fitted curve not only within the interval 3 < E < 4.5 covered by 
the "experimental" points, but also at the nearby points to the left and to the right of that 
interval. This is because we use proper analytic structure of the fitting S'-matrix. 

With thus found set of parameters, we located two zeros (0) of the function (ITUl) nearest to 
the real axis. They gave us the approximate (recovered from experimental data) resonance 
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energies: 



= 3.426388 - ^0.025531 , (20) 
£f = 4.821657-^2.036732. (21) 

Comparing them with the corresponding exact values (TTH!) and (1190 . we see that the result 
of fitting is very accurate for the S-matrix not only on the real axis but also in the nearby 
domain of the complex .E-surface. Actually, when choosing the "experimental" points we did 
not intend to reproduce the second resonance which is rather far away from the real axis. 



In addition to finding the optimal values of the parameters, the minimization program "Ml- 
NUIT" provides statistical errors (standard deviations) for them. This is done by calculating 
the matrix of partial derivatives of the minimized function with respect to all the parame- 
ters [231124]. Using these statistical errors, we can estimate the corresponding errors of the 
resonance parameters we found. 

To this end we considered the optimal values of the parameters as their mean values and 
randomly varied all the parameters around these values using a random-number generator 
with the Gaussian distribution of the width equal to the statistical errors. For each random 
choice of the parameters, we located zeros of the Jost function and then calculated their mean 
values and the standard deviations. After 1000 variations, we obtained: 

£f = (3.426309 ±0.005361)-^ (0.025921 ± 0.010444) , (22) 

= (3.810114 ± 0.448626) - - (1.457056 ± 1.090140) . (23) 

As one would expect, the reliability of the recovery of the second resonance is not that good 
as for the first one (around which the "experimental" points were taken). We however did not 
even expected to recover the second resonance at all. 



5.2 Two-channel model 

The two-channel potential. 



V{r) = ( 1 r^e-^ , (24) 
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of famous Noro and Taylor model [25] extends the single-channel potential of Sec. 15.11 It 
is written in the same dimensionless units with equal reduced masses fii = ^2 and angular 
momenta £i = £2 = in both channels. The threshold energies for the channels are Ei = 
and E2 = 0.1. 

The first three resonances of the Noro- Taylor model are given in Table [T] (the calculations can 
be found, for example, in Ref. [2J). The other resonances supported by the potential (!24l) . are 
too wide and therefore too far from the real axis. 

Since in the previous section we already demonstrated how the method works for a narrow 
resonance, considering the two-channel model, we focus our attention on the second resonance 
of the Table [H which is rather wide. We took "experimental" data in the elastic channels 
1 — >^ 1 and 2 — )■ 2 within the energy interval 5 < < 9. If one considers this segment of the 
real axis as the diameter of a circle in the complex plane, then such a circle will include the 
resonance point which we are looking for. The artificial data points (25 in the 1 — )■ 1 and 25 
in the 2 — )■ 2 channels) are shown in Figs. [3]and|4l 

In the same way as for the single-channel model, these 50 data points were fitted using the 
approximate matrices (11^ with M = 5 and real Eq = 7.25. In such a case the matrices (11^ 
are real and in total we have to adjust 48 parameters (matrix elements of a„ and bn) in order to 
minimize the generalized x^-function (ITBf) . After a thousand attempts with randomly chosen 
initial values for a„ and 6„, the best value of the minimum we found was X'^ = 1.9 x 10^"^. 

The results of this fitting are shown in Figs. [5] and [6] (thick curves). With rather small value 
of the X"^ achieved, there are no visible difference between the exact (thin curves) and fitted 
curves within the interval covered by the "experimental" points. However the extremely sharp 
(first) resonance to the left of this interval is missing since there are no data points reflecting it. 

Using the same approximate (fitted) Jost matrices, we calculated the cross sections for the 
transition processes 1 — )■ 2 and 2 — > 1, for which no data points were taken. In Figs. [7]and[Hl 
the comparison of the approximate (thick) and exact (thin) curves shows that we are able to 
rather accurately predict the cross section in one channel of the reaction on the basis of the 
data available in the other channels. 

Looking for the zero (nearest to Eq) of the determinant of the approximate Jost matrix, we 
found the following resonance: 

= 7.250742 - -1.513332 , (25) 
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which is very close to its exact location (see the second line of Table [T]). 



In order to find the partial widths, we use the method described in Ref. [2J, where it was 
shown that the ratio of the partial widhts can be found using the matrix elements of the Jost 
matrices, namely. 



p(out) p(in) _ p{out) p{in) 
-^11 ^22 ^12 ^21 



p{out) p(in) _ p{out) p(in) 
^22 ^11 ^21 ^12 



(26) 

E=S 



Together with the fact that Fi + r2 = F, the knowledge of the Jost matrices allows us to 
easily find Fi and F2. For the result ( 1^51) . this gives 

Ff = 0.347732 , F^* = 1.165600 . (27) 

These values reasonably well reproduce the corresponding exact partial widths given in Tabled] 

In the same way as for the single-channel case, we randomly varied (1000 times) the optimal 
parameters a„ and 6„, using the errors provided by the "MINUIT", and found the following 
mean values of the resonance energy and widths together with the corresponding standard 
deviations: 

= (7.251593 ± 0.292781) - ^ (1.085568 ± 0.537132) , (28) 

rf = 0.346348 ± 0.343585 , F^* = 0.739220 ± 0.516272 . (29) 

Although the third resonance is too wide and far away from the point Eq, we made an attempt 
to locate it using the same optimal expansion parameters. What we found, 

E^^ = 8.857247 - ^2.847466 , 

significantly underestimates the width. This is not surprizing, of course. The truncated 
series (IT31) can only be accurate within certain circle around Eq. In order to recover the third 
resonance one has to either take more terms in the series ( IT3l ) or shift the point Eq (the center 
of expansion) down in the complex plane. In both cases the number of fitting parameters would 
increase. 



6 Conclusion 

The proposed method is based on proper analytic structure of the parametrized S"- matrix that 
is used to fit experimental data. This mathematical correctness guarantees that after fitting 
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the data, we obtain the S'-matrix that is valid in all the channels, even in those where no data 
are available. This means that we can obtain the cross section for the channels, which are 
experimentally inaccessible, using the data in the other channels. 

The S-matrix properly fitted to the data for real energies is also valid at the nearby complex 
energies. This enables us to extract the resonance parameters as the real and imaginary parts 
of the zeros of the Jost matrix determinant that coincide with the S'-matrix poles. In addition 
to the total F, we are able to rather accurately obtain the channel partial widths. 

The most important limitation of the method described in this paper, is the fact that in its 
present form the method is only applicable to the systems with short-range interaction forces. 
A rigorous extension of the method that would include the Coulomb forces, could be done in 
a way similar to the one described in Ref [16J. This however would require a modified, much 
more complicated expression for the Jost matrix where all the non-analytic factors (square-root 
and logarithmic branching points etc.) are explicitly factorized. 

The other limitation is the non-relativistic character of the theory used to construct the S'- 
matrix. Although the method can still be used in a very wide range of problems dealing with 
low-energy atomic and molecular collisions, its possible applications in the intermediate- and 
high-energy particle physics would require relativistic corrections. This could be done in a way 
that is customary for those who work with mesons, namely, via using relativistic kinematics. 
For example, the non-relativistic channel momenta in the Jost matrices (ITUl [Tl]) can be 
replaced with the corresponding relativistic ones. 



If one accepts the approximate approaches that are traditional for meson-nuclear physics, then 
using such relativistic momenta in our S'-matrix that has the correct analytic structure, would 
be appropriate. Then in a similar simplified fashion the Coulomb corrections could also be 
introduced as it is done, for example, in Ref. [26J, where the /C-matrix for a short-range 
interaction is simply multiplied by a Coulomb-barrier factor Cq(?7) = 27X7] /{e^'^''' — 1). In our 
case, this trick is equivalent to multiplication of the matrices A and B in Eqs. (llOlllip by 
Co(r7) and l/Co(r7), respectively. The logic behind this is that for charged particles the Riccati- 
Bessel and Riccati-Neumann functions j£(A;r) and ye{kr) are replaced with the corresponding 
Coulomb functions Fi(r], kr) and Ge{r], kr), which at short distances differ from ji and ye by 




where -Ekm = c 



a/p^ -|- /U^c^ 



fic^ is the kinetic (collision) energy. 
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the factors and C'[ , and the matrices A and B are the factors combining and into 
the regular solution (see Eq. (41) of Ref. [16]). Of course this kind of relativistic and Coulomb 
corrections are not rigorously justified. However, their usage in meson-nuclear physics is based 
on some reasonable intuitive argumentation and proved to be working in practical applications. 
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Figure 1: Artificial data points for the single-channel model ( ITTI) . 
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Figure 2: Exact cross section (thin curve) for the single-channel model ( ITT!) and the result 
of fitting of the data (thick curve) with M = 5 and Eq = 3.4 in Eq. ( O . 
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Figure 3: Artificial data points for the first elastic channel of the model ( [24|) . 




Figure 4: Artificial data points for the second elastic channel of the model (!24l) . 
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Figure 5: Exact elastic cross section 1 — )• 1 (thin curve) for the two-channel model ( 1241) and 
the result of fitting of the data (thick curve) with M = 5 and Eq = 7.25 in Eq. ( IT3ll . 
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Figure 6: Exact elastic cross section 2—^2 (thin curve) for the two-channel model ([24]) and 
the result of fitting of the data (thick curve) with M = 5 and Eq = 7.25 in Eq. (11^ . 
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Figure 7: Exact inelastic cross section 1 — t- 2 (thin curve) for the two-channel model ([24]) 
and the prediction (thick curve) based on fitting of the data in the elastic channels 1 — )• 1 and 
2 -> 2 with M = 5 and Eo = 7.25 in Eq. (HH]). 
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Figure 8: Exact inelastic cross section 2 — )■ 1 (thin curve) for the two-channel model (^41) 
and the prediction (thick curve) based on fitting of the data in the elastic channels 1 — )■ 1 and 
2^2 with M = 5 and Eq = 7.25 in Eq. (HH]). 
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